Mesoscopic ensembles of polar bosons in triple-well potentials 
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Mesoscopic dipolar Bose gases in triple-well potentials offer a minimal system for the analysis of 
the non-local character of the dipolar interaction. We show that this non-local character may be 
clearly revealed by a variety of possible ground-state phases. In addition, an appropriate control of 
short-range and dipolar interactions may lead to novel scenarios for the dynamics of polar bosons 
in lattices, including the dynamical creation of mesoscopic quantum superpositions, which may be 
employed in the design of Heisenberg-limited atom interferometers. 
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Interparticle interactions are crucial in quantum 
gases pQ . They can usually be described by a short-range 
isotropic potential proportional to the scattering length 
a. Recently, dipolar quantum gases, in which the long- 
range and anisotropic dipole-dipole interaction (DDI) be- 
tween magnetic or electric dipole moments plays a signif- 
icant or even dominant role, have attracted a lot of inter- 
est as they show fascinating novel properties [2j [3: . To 
date, dipolar effects have been observed experimentally 
only with atomic magnetic dipoles, being particularly rel- 
evant in Bose-Einstein condensates (BECs) of 52 Cr where 
exciting new physics has been observed [U|Z] • Dipolar ef- 
fects have also been reported in spinor BECs [8], and in 
39 K and 7 Li BECs with a = HD]. Recent experi- 
ments with polar molecules [TTJ [12] open fascinating per- 
spectives towards the realization of highly-dipolar gases. 

Although a very clear and direct demonstration of the 
anisotropy of the DDI was given by the d-wave collapse of 
a Cr-BEC [6j[7], an equivalently obvious 'visual' proof of 
the non-local character of the DDI is still missing. Such 
a non-ambiguous qualitative evidence of the non-local 
character of the dipolar interaction could be provided in 
principle by the observation of novel quantum phases (su- 
persolid, checkerboard) in optical lattices [13] . However, 
the unambiguous detection of such phases is far from 
trivial, as is the preparation of the ground state of the 
system due to a large number of metastable states [T4] . 

In this Letter, we investigate a minimal system, namely 
a mesoscopic sample of dipolar bosons in a triple-well 
potential, which minimizes these restrictions, while still 
presenting clear visual non-local features (see "phase" B 
below). Non-dipolar BECs in double- well potentials have 
allowed for the observation of Josephson oscillations and 
non-linear self trapping [15], showing clearly that 'slic- 
ing' a BEC dramatically enhances the effects of interac- 
tions. The two- well Josephson physics is affected quanti- 
tatively (although not qualitatively) by the DDI [16] [17] 
(the DDI may induce however significant inter-site ef- 
fects in coupled ID and 2D bilayer systems [18-20 ). On 



the contrary, as we show below, the DDI does introduce 
qualitatively novel physics in the Josephson-like dynam- 
ics in three- well systems. We discuss how the DDI leads 
to various possible ground states, which may visually re- 
veal the non-locality of the DDI. In addition, we show 
how this non-locality leads to a peculiar quantum dynam- 
ics characterized by striking new phenomena, including 
the dynamical formation of mesoscopic quantum super- 
positions (MQS). MQSs produced in cavity QED or with 
trapped ions [21] require complex manipulations, whereas 
in the present system, they arise naturally, similar to the 
MQSs obtained in BECs with attractive interactions in 
double wells [22, 23] or lattices [24]. We then comment on 
the design of four-site Heisenberg-limited atom interfer- 
ometers using the dynamical creation of MQS, and finally 
discuss possible experimental scenarios. 

We consider N dipolar bosons in a three- well potential 
Vtrap(^) (Fig. [T^l). The wells are aligned along the ?/-axis, 
separated by a distance £ and an energy barrier Vq. The 
bosons are polarized by a sufficiently large external field, 
with a dipole moment d along a given direction. The 
lattice potential is strong enough compared to other en- 
ergies (in particular the interaction energies) such that 
the on-site wavefunctions (j)i=i^^{ r ) are fixed, being in- 
dependent of the number of atoms per site. For a large- 
enough Vo we may assume 4>i(r) = (j)(r — r$), where Vi is 
the center of site i. In addition we may assume <f> to be a 
Gaussian with widths cr X} y yZ . We limit to the case where 
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FIG. 1: (color online) (a) Schematic view of the three- well 
system; (b) MQS interferometer with four wells (see text). 
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a y is small enough with respect to t so that the sites are 
well defined. Re-expressing the bosonic field operator as 
VK r ) = Y^h=i we ma y write the Hamiltonian as: 



H 



-J 



a\ (ai + a 3 ) + h.c. +^^7^(7^-1) 



h\h2 + h 2 h 3 H — nin 3 



(1) 



where J = - / dr^r) [-h 2 V 2 /2m + Vt rap (r)] ^ 2 (r) 
is the hopping rate, Uq = g J|0i| 4 dr + 
/ |0i(^)| 2 | ( / ) i(^ / )| 2 ^ r dd(^ — r')&r&r' characterizes 
the on-site interactions, U\ = J \(j)i(r)\ 2 \(j) 2 (r f )\ 2 Udd(^ — 
r')drdr' is the coupling constant for nearest-neighbor 
DDI, and hj = OjCtj. In the previous expressions 
g = 47rh 2 a/m is the coupling constant for the short- 
range interactions, with a the s-wave scattering length. 
The DDI is given by Udd(r) = d 2 (l — 3 cos 2 #)/r 3 , where 
6 is the angle between r and d, d 2 = fiofi 2 /(47r) for 
magnetic dipoles (/i is the magnetic dipole moment) or 
d 2 = J 2 /47reo for electric dipoles (J is the electric dipole 
moment). The parameter a in Eq. depends on the 
geometry of Vtrap(V) (a = 8 if the wavefunctions are well 
localized in all directions compared to i, and decreases 
towards a = 4 when a x /£ --> 00 [25 J. In the following 
we focus on the localized case, i.e. a = 8, but all results 
remain valid for 4 < a < 8. Finally, note that Uo results 
from short-range interactions and DDI, and that the 
ratio between Uo and U\ may be easily manipulated by 
means of Feshbach resonances, by modifying the dipole 
orientation d, and by changing t [25] . 

Since J2i hi = N is conserved by ([!]), we may re- write 
H (up to a global energy UoN(N — l)/2) as an effective 
Hamiltonian without on-site interactions: 



H 



a\ (ai + a 3 ) + h.c. + (Ui - U )n 2 [fi\ + h 3 ] 
- Uo) funs. (2) 



The gross structure of the ground-state diagram is un- 
derstood from the J = case, where the Fock states 
|ni, 77,3) are eigenstates of H, with energy £?(ni, n 3 ) (since 
N is conserved, the Fock states are defined by ni >3 ). 
The minimization of E provides four classical "phases". 
For U > and U x < 8/7 /15, and U < and 
U\ < —8\Uo\, phase (A) occurs, with 77-1 = n 3 = [n/2\ 
with n = 16N(U - U 1 )/(24U - 31Z7i) (where |_-J de- 
notes the integer part). Phase (B) appears for Uo > 
and 8Z7o/15 < U\ < 8/7o, being characterized by n\ = 
n 3 = 7V/2. For C7 > and Ui > 8U 0: and U < and 
/7i > —\Uo\ phase (C) occurs, with n 2 = N (actually 
states with m = N are degenerated, but the degeneracy 
is broken by tunneling which favors n 2 = N). Finally, 
phase (D) occurs for Uq < and 8U0 < U\ < Uq, being 
characterized by a broken symmetry, with two degener- 
ated states with n\ = [n\ , n 3 = and vice- versa. 



Fig. [8ja) shows (n)/N, with ft = hi + n 3 for N = 18. 
We can see that phases (A)-(D) describe well the gross 
structure of the ground-state diagram (a similar graph 
shows, as expected, that the (D) phase shows large fluctu- 
ations AS in S = h\ — n 3 ). However, tunneling is relevant 
at low \Uo\ and \U\\ and at the phase boundaries. In gen- 
eral, the system is in a quantum superposition of differ- 
ent Fock states |^) = $X 1=0 ^n^o C ( n ii n s)\ n i, n s). 
Fig. ^b) depicts AS in the region C7o,i > 0. As expected 
at small \Uo,i\/J tunneling dominates and the product 
state (aJ/V^ + ^J/^ + a 3 /v // 2) Ar |vac) is retrieved (|vac) is 
the vacuum state). This state transforms into phase (A), 
which for growing Uo becomes the Fock state |iV/3, N/3). 
Phase (C) remains the Fock state |0, 0) (n 2 — A/"), and the 
border (B)-(C) is characterized by a first-order "phase 
transition" [26], at which n 2 abruptly jumps from to 
N. Fig. [8jc) represents schematically phases (A) to (D). 

Phase (B) is characterized by vanishing (h 2 ) and An 2 , 
and (hi) = (n 3 ). It strikingly reveals the non-local char- 
acter of the DDI, similarly to the biconcave BECs pre- 
dicted in [27 , but with a much higher "contrast". Note 
however that the actual ground-state may significantly 
depart from | N/2, N/2), since |A<5| is significant at the 
(B)-(C) transition (Fig. gb)). At U x = 8U , the ground 
state is a coherent state (a\ + a^^lvac), i.e. coherence 
between the two extremal sites is preserved in spite of 
the absence of particles in site 2. This coherence is un- 
derstood from (11), since for Ui = 8U0 there is no effec- 
tive interaction between sites 1 and 3. Since (h 2 ) <C 1 
due to the effective repulsive nearest-neighbor interac- 
tions (Ui — Uo > 0), then sites 1 and 3 form an effective 
non-interacting two-well system coherently coupled by a 
second-order process through site 2 (with effective hop- 
ping J e ff = J 2 /7(N — l)Uo). Hence the coherent region 
extends inside (B) for \Ui — 8Uo\ < J e ff- Thus for larger 
NUo the coherent region shrinks (reducing to the very 
vicinity of Ui = 8U0 as seen in Fig. [8^b)). 

Such a 1-3 coherence has important consequences for 
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FIG. 2: (color online) (a) (h)/N as a function of £/o,i for 
N = 18; (b) AS/N in logarithmic scale for Uo,i. The dashed 
lines show the boundaries between the classical phases (A)- 
(D), that are shown schematically in (c). 
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the quantum dynamics, best illustrated by consider- 
ing initially all particles at site 3. Interestingly, (711,3) 
show perfect Josephson-like oscillations (with frequency 
2 J e ff /K) although for any time (77,2) = An2 <C 1. How- 
ever, J e ff decreases with NUo and hence the observation 
of this effect demands a mesoscopic sample, since other- 
wise the dynamics may become prohibitively slow. Off 
the U\ = 8/7o boundary, inside phase (B), the residual 
1-3 interaction leads to a damping of the Josephson os- 
cillations (connected to number squeezing). Eventually 
for \U\ — 8t/o | ^> J e ff self-trapping in 3 occurs. 

Phase (D) is characterized by a large AS and (77,2) 7^ 0, 
and two degenerated states: = (i), and n\ = 
(ii). Strictly speaking the exact ground state is provided 
by a MQS of these two states, but the gap between the 
ground-state and the first excited one is vanishingly small 
(<C J) even at the U\ = Uo < boundary and for N 
as small as 18. Experimentally, the signature of phase 
(D) would thus consist in measuring large shot-to-shot 
fluctuations in <5, while never observing simultaneously 
atoms in both sites 1 and 3. At U\ = Uo < 0, states 
(i) and (ii) become coherent superpositions of the form 
(a\ + a\) N I vac) and (a\ +ag) iV |vac), respectively. These 



superpositions may be understood from Eq. (11), which 
for U\ = Uq < becomes 



H = — J a\ (ai + as) + h.c. 



7|tfo I 



-nin 3 . 



(3) 



which describes a non-interacting two-well system if n\ = 
or ris = 0, leading to the coherent states (i) and (ii). 

Hamiltonian (|3| leads to an intriguing quantum dy- 
namics characterized by the creation of MQSs. From an 
initial Fock state |0,0) (ri2 — A"), if a particle tunnels 
into site 1 (state |1,0)), a subsequent tunneling from 2 
to 3 (state |1,1)) is produced with a bosonic-enhanced 
hopping rate J\/N — 1. However, the state |1, 1) has an 
interaction energy 7|f7 |/8. Hence if J <C 7\Uq\/%\/N - 1 
then the tunneling from 2 to 3 remains precluded. On the 
contrary the hopping into 1 presents no energy penalty. 
As a result, if the first particle tunnels into 1, then a co- 
herent 1 — 2 superposition is established. Of course if the 
first particle tunnels into 3, then a 2 — 3 superposition 
occurs. Since the initial process is coherently produced 
in both directions, then a MQS |$(t))|0) + |0)|$(t)) is 
formed, where \<&(t)) = Yln=o C n (t)\n), with the normal- 
ization condition 2^ =1 |C n (£)| 2 + 4|C (t)| 2 = 1 [25]. 
Fig. [3^1 shows that (hi^)(t) perform a coherent oscil- 
lation, which however damps for longer times. This 
damping is again a remarkable consequence of the non- 
local character of the DDI. Virtual hoppings of a sin- 
gle particle from site 2 into site 3 (1) induce a second- 
order correction of the energy of the states |n, 0) (|0, n)): 
AE n = SJ 2 (N - n)/7\Uo\n, which distorts the Joseph- 
son Hamiltonian, and leads to a significant damping after 
a time scale of the order of r ~ 7|£7o|/8J 2 A" (in agree- 
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FIG. 3: (color online) (a) (ni, 3 (t)) (dashed), (n 2 (t)) (solid), 
for Uo = Ui = -100 J and N = 18; (b) Probability P 4 (A) as 
a function of (j) for the interferometric 4-site arrangement (see 
text) with N 14, U = Ui = -100 and Jt 2.7. 



ment with our numerics) [25 . At longer times, chaotic 
dynamics may even occur [28] . 

The three- well system hence acts as a MQS-splitter un- 
der the mentioned conditions. We stress, however, that 
a MQS (although asymmetric) is still created [25 , even 
for unequal hoppings Jij for nearest neighbors, as long as 
Ji2,23 < 7\U \/8y/N -1. We note also that if U x ^ U 
a MQS is created if \U\ — Uo\ ^ J, but nearest-neighbor 
interactions enhance the damping in each MQS branch. 
If \U\ — Uo\ ^> J bosons at site 2 remain self-trapped. 

The MQS-splitter opens fascinating possibilities be- 
yond the 3-well system, most relevantly in the context 
of Heisenberg-limited atom interferometry. We illustrate 
this possibility by considering a simple interferometer 
based on a four- well system (Fig.JlJ)). Initially the bosons 
are at site 2 (which acts as the input port). Sites 1 and 

3 play the role of the interferometer arms, whereas site 

4 acts as the output port, where the interferometric sig- 
nal is read out. We consider hoppings J21 = J23 = J, 
but J34 = Je 1 ^ = J^ 4 . We are interested in the <\>- 
sensitivity of the population at site 4. This arrangement 
is chosen for its theoretical simplicity (more general ar- 
rangements work along similar lines), although it may 
be implemented also in practice by means of Raman- 
tunneling [29 . Under the MQS conditions (in this case 

l)|C/ |/2x/2), the 



U x = Uo < and Jv^A^T <C (2y/2 
system evolves into an entangled MQS formed by Fock 
states such that riifij = for next-nearest neighbors. It 
is straightforward to show that the probability to find 
Af particles at site 4 depends explicitly on the phase <\> 
as P 4 (A0 ~ cos 2 (A'0) (P 4 (n ^ N) are only indirectly 
^-dependent due to normalization). Hence P^N) has a 
modulation of period 5(f) = ir/N (Fig. [3]o) , contrary to 
the period 5(j) = tt expected for independent single par- 
ticles, allowing for a Heisenberg-limited interferometric 
measurement of the phase <fi. This super-resolution is an 
unambiguous signature of the coherent character of the 
MQS thus created [30j|3j. (77,4) presents a similar mod- 
ulation (but with poorer contrast). Calculations with a 
six-site arrangement provide similar results [25] . 

In the final part of this Letter we discuss experimen- 
tal feasibility. Triple- well potentials as in Fig. [l] may be 
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controllably implemented with optical potentials. By su- 
perimposing, onto a single-beam optical trap which pro- 
vides the ^-confinement, a tightly focused beam (with 
a waist ~ 1/im, see e.g. [32 J, one may create a tight 
'dimple' acting as one well. To realize a triple-well (or 
even more complex configurations), several possibilities 
exist. Using an acousto-optic modulator (AOM) with 
several rf frequencies [33j [34], several diffracted beams 
are created, whose intensity and position can be con- 
trolled independently. Another option using an AOM 
consists in toggling the dimple between several positions 
at high rate, to create almost arbitrary time- averaged 
potentials [35]. Such an implementation has several ad- 
vantages: arbitrary, time-dependent energy offsets can 
be applied to the different sites; the inter-site separation 
£ can be changed in real time, easing the preparation of 
a given atom number in each well (e.g. by performing 
evaporative cooling with different energy offsets in each 
site), and the detection of the population in each well 
(before imaging, Vq may be increased to freeze out the 
dynamics and then £ increased, thus relaxing constraints 
on the imaging resolution). 

We now evaluate J, Uo and U\ for realistic experimen- 
tal values. Although in our calculations we have just con- 
sidered N up to 36, similar ground- states are expected 
for larger N (but, as mentioned above, the observation of 
the quantum features at the (B)-(C) and (D)-(C) bound- 
aries demands small samples). In particular, consider a 
triple-well potential formed by three Gaussian beams of 
waist 1 fim separated by £ = 1.7 /mi. For a barrier 
height Vo/h ~ 2500 Hz, we obtain J/h ~ 10 Hz, and 
the typical value of NU ± /J is then - 10 for N = 2000 
52 Cr atoms. The value of Uo can be tuned, for a fixed 
geometry, by means of Feshbach resonances [4], so that 
one can explore e.g. the first-order (B)-(C) "transition" 
with 52 Cr by varying U\/Uq. However, the MQS cre- 
ation demands small samples, being hence more realistic 
with polar molecules. For example, for KRb molecules 
placed at a distance £ = 1 /im and maximally polarized 
(d = 0.5 D) parallel to the joining line between the sites, 
Ui/h ~ — 70Hz. Under these conditions the MQS con- 
dition implies, for N = 36 molecules, J/h of a few Hz. 
Single-atom sensitivity has been achieved with fluores- 
cence imaging [36], so that the relatively small values of 
N considered here should be detectable. 

In summary, we have studied a simple system of dipo- 
lar bosons in a triple well, showing that the non-locality 
of the DDI leads to qualitatively novel physics that may 
be explored with a high degree of control over all pa- 
rameters via the trap geometry, dipole orientation, and 
Feshbach resonances. We have shown that the ground- 
state phases present abrupt crossovers induced by the 
non-local nature of the DDI, which may be explored with 
52 Cr BECs. In addition, the dynamics presents intriguing 
new scenarios, especially for the case of polar molecules, 
including the dynamical creation of MQSs, which may be 



employed for Heisenberg-limited interferometry. 
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Supplementary material for 
Mesoscopic ensembles of polar bosons in triple-well potentials" 



In this supplementary text, we provide additional de- 
tails concerning several points discussed in the Letter. 



VARIATION OF THE PARAMETER a WITH 
THE RATIO £/a x 



The parameter a is defined as the ratio 



a 



E(a x ,2£) 



(4) 



where E(a x ,£) is the dipolar interaction energy of two 
Gaussian clouds of length a x along the x-axis, containing 
N particles, and separated by a distance £ along the y 
direction. Here we consider that the dipoles are pointing 
towards the z direction as in Fig.[4^a) (implying U\ > 0), 
but the results hold for other dipole orientations. We 
assume that the sizes a VyZ of the cloud along the y and z 
axes are small as compared to I. Therefore, E depends 
only on i and a x . From scale invariance, we deduce that 
a depends only on the ratio £/<j x . We can easily work out 
two limits: (i) if a x <C £, the two clouds can be considered 
as point-like, and the dipolar energy of the two clouds 
thus scales as E ~ l/£ 3 . In that case, one thus gets, 
from Q, a = 8; (ii) in the opposite case a x ^> £, one 
practically deals with two infinite chains of dipoles, and 
the interaction energy is then known to scale as E ~ l/£ 2 , 



t t t 

(a) U > and Ui > 

z 

)—*y 



.4, 



4 < 4 

(b) U < and U x > 



(c) U > and U Y < 



(d) U o <0 and U v < 



FIG. 4: Changing the geometry of the sites and the orienta- 
tion of the dipoles allows, for three aligned wells, to choose 
any combination of signs for Uq,Ui. For cases (b) and (c) , the 
signs given on the figure are correct only for moderate aspect 
ratios of the on-site wavefunctions. 



yielding a = 4. In between these two cases, a varies 
monotonously with the ratio £/cr x , as we shall see below. 

Since we assume a y ^ z small compared to £, we can ap- 
proximate the Gaussian atomic densities by: 



rii{x,y,z) = Aexp 



ri2(x,y,z) = Aexp 



S(y)S(z), 



5(y-£)S(z), 



where A is a normalization constant. Calculating the 
dipolar interaction energy 

E = J dr dr' ni{r)n2{r , )U^{r — r')^ 



we get 



E(a x ,£) ex dxdx 



, exp [— (x 2 



\(x - x') 2 + £ 2 ] 



3/2 



(5) 



From this expression, it is clear that a depends only on 
the ratio £/a x , as expected from scale invariance; this 
dependence is shown on Fig. [5j where the curve was ob- 
tained by a numerical integration of ([5|. 



TUNING OF THE INTERACTION 
PARAMETERS U AND Ui 

In this section we elaborate on how the interaction 
parameters can be modified such that the sign of Uo,i 
can be tuned. First, let us consider the arrangement of 
three aligned wells. Disregarding the contact interaction 
contribution to E/o? all sign combinations for C7o,i can be 
achieved by exploiting the anisotropy of the DDI in a 
carefully designed geometry, as sketched in Figs. [4j 




FIG. 5: Parameter a as a function of £/a x , obtained from |5]). 
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We now turn to the the square (or even hexagonal) 
arrangement for the interferometric setup (see the Letter 
and sections below), for which we want to have Uq < 
and XJ\ < 0. In this case, because of geometric con- 
straints, the sign of the DDI must be inverted. In prin- 
ciple, this can be achieved for magnetic dipoles by using 
a rotating field pQ (although this precludes the use of 
magnetically induced Feshbach resonances). In the case 
of polar molecules (more relevant for the study of the 
MQS splitter as discussed in the Letter), dressing the 
rotational levels by AC fields, as proposed in [2], allows 
for inverting the sign of the DDI; in combination with a 
proper design of the trapping geometry, this should allow 
for the realization of the proposed experiments. 
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DYNAMICS OF THE MQS SPLITTER 



In this section we comment on the dynamics of the 
MQS-splitter in detail. Under the MQS conditions U\ = 
U < and Jy/N - 1 <C 7\U \/S, only states of the form 
(using the notation of the Letter) |n, 0) or |0, n) are pos- 
sible. Under these conditions the system (originally pre- 
pared in 1 0, 0) ) evolves dynamically into a MQS state of 
the general form: 



Wt)> = |$>|0) + |0)|*), 



(6) 



with |$> = £„ =0 C„(t)|iV - n), where the C n {t) 
coefficients fulfill the normalization condition 

2En=0,N-l\Cn(t)\ 2 + MC N (t)\ 2 = l. 

Note that a slightly different tunneling J 12 = J + S J, 
J23 = J-SJ such that both J 12 ^VN - 1 <C 7|U |/8 will 
lead to a similar state, although now slightly asymmetric: 



|V(t)> = |*>|0) + |0>|*'>. 



(7) 



In Fig. [6] we consider the case J12 = 0.9 J and J23 = 1.1 J, 
Uq = U\ = —100 J. One may see that the dynamics 
is, as expected, different for sites 1 and 3 (this should 
be compared with Fig. 3(a) in the Letter). However, 
the state is still of the MQS form, since (hihs) remains 
negligible (< 3 x 10 -3 ) at any time. 

Let us go back to the symmetric case J12 = J23 = J- 
Fig. 3b in the Letter shows the evolution of the pop- 
ulations for Uo = —100 J = J7i, and N = 18 particles. 
We have evaluated that the MQS condition is actually 
fulfilled, i.e. (77472,3) remains negligible at all times. The 
population in 1 follows a coherent Josephson oscillation, 
but it damps after several oscillations (for Uo = —10 J, 
only one oscillation remains). The reason behind this 
damping is interesting, since it constitutes itself a quite 
remarkable non-local effect introduced by the dipole- 
dipole interactions. 

The main reason behind this damping is the existence 
of second-order processes, coming from virtual transi- 
tions of the form |0,n) — » |l,n) — >• |0,n), i.e. a single 



particle tunnels to the left from site 2 and comes back. 
These virtual explorations of the left side (by a single par- 
ticle!) result in an energy shift for the state |0,n), which 
acquires a second order energy shift 8 J 2 (TV — n)/7\Uo\n. 
This means that there is a coupling between different 
many-body states which have different energies. The 
physics of the 2-3 system is then given by the effective 
Hamiltonian (here \n) means only the site 3) 



N 

^eff ^ -J^yfti[\n-l)(n\ + \n)(n-l[ 



N 



8J 2 (N-n) 



n)(n\. 



(8) 



Note that without the extra second-order correction we 
have the usual Josephson coupling (no interactions). But 
with the second-order term the problem is now truly 
many-body, and we cannot write the Hamiltonian any 
more as a Josephson-like Hamiltonian. This second-order 
shift is the responsible of the observed damping. This 
may be easily seen by having a look to the typical en- 
ergy scale of the perturbation ~ 8J 2 7V/7|C/o| — 0.18 J 
(for the example of Fig. 3b of the Letter), which leads to 
a damping time scale of the order of 6/ J, which is what 
one observes in the numerics. 

Summarizing, the system behaves indeed as a MQS- 
splitter. Second-order processes (due to single-particle 
tunnelings to the left for right moving states, or viceversa 
for left-moving states) lead to a damping of the Joseph- 
son oscillations, and hence induce a complex and rich 
dynamics for the MQSs. This has important practical 
consequences, e.g. there are optimum times for the in- 
terferometric measurements (see below). 
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INTERFEROMETRIC ARRANGEMENT USING 
A FOUR-SITE SQUARE SYSTEM 

In this section we would like to discuss the use of the 
MQS-splitter idea in the context of Heisenberg-limited 
interferometry. 

We consider a 4-site arrangement as in Fig. 1(b) of 
the Letter. There is tunneling between nearest neighbors 
with hopping constant J. However, we shall introduce an 
additional phase </>, such that the tunneling J 3 ^ 4 — e l( ?J 
and Ji^4 = e~ l( ?J. We are interested in the sensitivity 
of the system to this phase. We will consider that the 
initial population is all in site 2. We shall now evaluate 
the ^-sensitivity of the population at site 4. 

The Hamiltonian of the system is of the form: H = 
Ht + Hi, where the hopping part is given by 



H T = -J 



a\(a\ + as) + h.c. 



j[4( ei ^3 + e _i %)+h.c.l , 



(9) 



and the interaction part is 
Hi = ^yJ2hi(hi-l) 



{hi + n 3 )(n 2 + n 4 ) + ^-^=(nin 3 + n 2 n 4 ) 



10) 



As for the triangle, we may eliminate a constant and get 



the effective interaction part: 

Hi = (Ui - Uo)(ni + n 3 )(n 2 + n 4 ) 



(11) 



Note that now the sites are at the vertices of a square, and 
hence the next-nearest-neighbor interaction is Ui/2y/2 
(as for the 3-site case we assume here point-like sites, see 
the discussion of the first section above). 

We may again find the MQS-splitter condition U\ = 
Uq < 0. In that case 



Hi 



2y/2-l 
2V2 



\Uo\(hin 3 + n 2 n 4 ). (12) 



We have again a very similar MQS-splitter condition as 
that in the Letter: 



'2y/2- 



2V2 



\Uo\- 



(13) 



When this is fulfilled, then starting from all the popula- 
tion in site 2, if 1 is populated 3 is not, and viceversa. 
Moreover, due to the (effective) repulsion between sites 
2 and 4, 4 may be only populated once 2 is empty. 

This means that the state of the system is of the form 
(with the Fock-state notation |ni, n 2 , n 3 , n 4 )): 



\1>{t)) = c(t)^D n (t){\N - n,nA0) + \0,n,N - n,0)) 

n 

+ s(t) J2 Cn(4>, t)(e in4, \0, 0,N-n,n)+ e~ in *\N - n, 0, 0, n», (14) 

n 

I 



with the normalizations: 

|c(t)| 2 + | S (*)| 2 = 1, (15) 
2^ \D n (t)\ 2 +4\D N (t)\ 2 = 1, (16) 

2^\Cn(^t)\ 2 ^4\C N (^t)\ 2 cos 2 N^ = 1. (17) 

In principle the <p dependence appears explicitly only 
for the n 4 = N term, but note that due to the normal- 
ization condition it appears also in terms with n 4 7^ N . 
The probability to find all the atoms at site 4 is hence 
P 4 (7V) = 4|Cjv((/>, t)\ 2 cos 2 Ncj), and hence it has a mod- 
ulation of period 5(j) = ir/N. The modulation of P 4 (iV) 
has clearly a contrast 1 (see Fig. 3(b) of the Letter). 
The population in 4 is hence very sensitive to the phase, 



and, if monitored, allows for a Heisenberg-limited inter- 
ferometric measurement of the phase <j). Obviously, if, 
instead of the above coherent superposition, one creates 
a statistical mixture, the probability to find all the par- 
ticles at site 4 is independent of <p. 

In Fig. 3(b) in the Letter we have chosen the time 
Jt = 2.7 for the interferometric measurement. This is 
done so, since for this time the transfer into site 4 is 
optimal, i.e. the probability Pa(N) is maximal for <fi = 
(see Fig. [7]). 

One may also monitor the average population (n 4 ). 
This population is also modulated (see Fig. [8| but the 
contrast is poorer. The reason may be found in the nor- 
malization condition, since when the population P 4 (iV) 
tends to zero, this increases the relative probability to oc- 
cupy n 4 7^ N. The latter may be observed by monitoring 
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FIG. 7: P±(N) as a function of Jt for N = 14 and Uo = Ui = 
-100J. 
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FIG. 8: Average population (fiA)/N as a function of Ntp/ir 
for the 4-site system with N = 14, Uo — Ui — —100 J after 
Jt — 2.7 (dashed red curve), and for the 6-site system (solid 
blue curve), at time Jt = 3.48 and same E/o,i- 



HEXAGONAL ARRANGEMENT 

In this last section we briefly discuss an hexagonal in- 
terferometric arrangement, showing that the Heisenberg- 
limited interferometry based in the MQS-splitter may be 
extended to more general setups. We consider the ar- 
rangement shown in Fig. [9] The full quantum simulation 
of the dynamics becomes much harder, and may be only 
performed for a rather small number of particles. We 
have performed simulations for up to N = 6 particles. 
Initially all the population is at site 1 (input port of the 
interferometer). Sites 2 — 3 and 6 — 5 act as the arms of 
the interferometer, whereas site 4 acts as the output port, 
where the interferometric measurement is performed. As 
for the 4-site arrangement all hoppings are the same, ex- 
cept for the hoppings J34 and J54 which have an extra 
phase. Pa(N) presents a perfect modulation with period 
5<j) = tt/N, very similar as for the 4-site system (Fig. 3(b) 
in the Letter). In Fig. [8] we show the comparison for 4 
and 6 sites of (n±)/N a a function of N(j)/7r. 




FIG. 9: Sketch of the 6-site hexagonal system employed as a 
simple model of an atom interferometer. 



the population distribution P±(n). Note also that, again 
due to normalization, {fi\^) are also modulated (out of 
phase) with period 5<fi = tt/N. 

These results must be compared with the non- 
interacting case Uo = U\ = 0, where Pa(N) shows a 
(trivial) modulation with period 5<j) = tt. 
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